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ABSTRACT 

We  present  observations  of  newly  discovered  24 /xm  circumstellar  structures  detected  with  MIPS  around  three 
evolved  stars  in  the  Cygnus-X  star-forming  region.  One  of  the  objects,  BD+43  3710,  has  a  bipolar  nebula, 
possibly  due  to  an  outllow  or  a  torus  of  material.  A  second,  HBHA  4202-22,  a  Wolf-Rayet  candidate,  shows 
a  circular  shell  of  24 /xm  emission  suggestive  of  either  a  limb-brightened  shell  or  disk  seen  face-on.  No  diffuse 
emission  was  detected  around  either  of  these  two  objects  in  the  Spiizer  3.6-8  /xm  IRAC  bands.  The  third  object 
is  the  luminous  blue  variable  candidate  G79.29+0.46.  We  resolved  the  previously  known  inner  ring  in  all  four 
IRAC  bands.  The  24 /xm  emission  from  the  inner  ring  extends  ^\'.2  beyond  the  shorter  wavelength  emission, 
well  beyond  what  can  be  attributed  to  the  difference  in  resolutions  between  MIPS  and  IRAC.  Additionally,  we 
have  discovered  an  outer  ring  of  24 /xm  emission,  possibly  due  to  an  earlier  episode  of  mass  loss.  For  the  two 
shell  stars,  we  present  the  results  of  radiative  transfer  models,  constraining  the  stellar  and  dust  shell  parameters. 
The  shells  are  composed  of  amorphous  carbon  grains,  plus  polycyclic  aromatic  hydrocarbons  in  the  case  of 
G79.29+0.46.  Both  G79.29+0.46  and  HBHA  4202-22  lie  behind  the  main  Cygnus-X  cloud.  Although  G79.29+0.46 
simply  may  be  on  the  far  side  of  the  cloud,  HBHA  4202-22  is  unrelated  to  the  Cygnus-X  star  formation  region. 

Key  words:  stars:  AGB  and  post-AGB  -  stars:  individual  (BD-»43  3710,  HBHA  4202-22,  G79.29— 0.46) 


1.  INTRODUCTION 

Cygnus-X  is  arguably  the  biggest  and  brightest  massive  star 
formation  region  within  2  kpc  of  the  Sun  (d  ~  1.7  kpc; 
Schneider  et  al.  2006,  2007).  It  contains  several  hundred  Hii 
regions,  numerous  OB  associations  with  over  a  thousand  OB 
stars,  and  thousands  of  low-mass  stars.  We  have  performed  an 
unbiased  survey  of  ~24  deg^  in  Cygnus-X  with  IRAC  (Fazio 
et  al.  2004)  and  MIPS  (Rieke  et  al.  2004)  on  the  Spiizer  Space 
Telescope  (Werner  et  al.  2004).  The  primary  objectives  of  this 
Spiizer  Legacy  project  were  to  study  the  evolution  of  high- 
mass  protostars,  clustering  of  both  high-  and  low-mass  stars  in 
the  complex,  as  well  as  several  other  goals  described  by  Hora 
et  al.  (2008b).  Follow-up  spectroscopy  with  Spiizer*^  Infrared 
Spectrograph  (IRS;  Houck  et  al.  2004)  was  obtained  for  several 
dozen  promising  sources  in  the  MIPS  and  IRAC  survey  regions. 

In  addition  to  the  numerous  young  stars  that  are  the  primary 
focus  of  the  Legacy  project,  there  are  dozens  of  evolved  objects 
in  the  survey  region,  from  Miras  and  Cepheids  to  Wolf-Rayet 
( WR)  and  carbon  stars  to  planetary  nebulae  (PNe)  and  supernova 
remnants  (SNRs).  After  the  initial  discovery  of  a  bipolar  nebula 


at  24  /xm  around  a  carbon  star  during  the  data  processing  quality 
control  effort,  a  systematic  search  for  additional  circumstellar 
nebulae  was  made.  SIMBAD  was  used  to  compile  positional 
lists  of  evolved  objects  that  lie  within  the  bounds  of  the  24 /xm 
observations.  The  24 /xm  images  were  then  inspected  for  the 
presence  of  point  sources  and  associated  nebulosity,  if  any. 
More  than  80%  (146)  of  the  191  known  evolved  objects  are 
detected  at  24  /xm  as  point  sources  (Table  1 ).  These  types  of 
objects  can  be  bright  in  the  infrared  for  a  variety  of  reasons  such 
as  emission  from  dust  or  forbidden  lines.  This,  of  course,  is  an 
incomplete  survey  of  the  evolved  stars  as  it  does  not  include 
the  faint  population  of  asymptotic  giant  branch  (AGB)  stars 
detected  with  Spiizer  (also  known  as  contaminants  in  the  young 
stellar  object  studies).  Be  that  as  it  may,  only  three  evolved 
.stars  in  the  Cygnus-X  region  were  distinguished  from  the  vast 
majority  by  exhibiting  extended  24  /xm  emission.  The  emission 
around  two  of  these  stars,  which  are  listed  in  Tabic  2,  was 
first  reported  by  Kraemer  et  al.  (2009).  Here,  we  discuss  the 
morphological  characteristics  of  the  circumstellar  structures  as 
well  as  the  spectral  properties  of  the  dust  and  the  central  exciting 
objects. 
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Table  1 

Evolved  Objects  in  Cygnus-X:  Point-source  Detections  at  24 


Object 

Type 

No.  in  24 /im 
Region 

Detected  at  24  /xm 
Region 

Detection 

Rate 

PNe 

7 

5 

71% 

WRs 

9 

9 

100% 

Post-AGBs 

2 

2 

100% 

Carbon  stars 

52 

51 

98% 

S  stars 

6 

5 

83% 

SNRs 

5? 

0? 

0?% 

Miras 

4 

4 

100% 

Semi-regulars 

4 

4 

100% 

Cepheids 

4 

4 

100% 

Dwarf  novae 

3 

1 

33% 

Other  V*s 

85 

61 

72% 

Notes.  Object  type  as  reported  in  Simbad.  This  is  likely  to  be  an  incomplete 
listing,  since  not  all  objects  are  typed  correctly.  BD+43  3710,  for  instance, 
is  listed  merely  as  a  star.  There  were  an  additional  40  variable  stars,  mostly 
generic  “V*”  type  but  including  one  RR  Lyr,  which  had  insufficiently  precise 
coordinates  to  determine  if  they  had  24  /im  counterparts.  It  is  not  clear  how  big 
the  SNRs  are  expected  to  be  at  24  /im,  so  the  data  inspection  may  not  have  been 
sensitive  to  their  emission. 

2.  OBSERVATIONS  AND  DATA  REDUCTION 

The  Cygnus-X  Legacy  Survey  (Hora  et  al.  2008b)  observed 
a  X  T  region  centered  on  {I,  b)  ^  (80°,  0^3)  with  Spitzer's 
IRAC  (all  four  bands)  and  MIPS  (24  and  70 /xm).  The  MIPS 
observations  followed  the  MIPSGAL  survey  strategy  (Carey 
et  al.  2009),  using  the  fast-scanning  mode  with  3  s  integrations 
and  148"  between  legs.  This  results  in  '^lO  times  redundancy 
per  pixel  or  ~30  s  total  integration  time.  We  use  the  24  pm  data 
processing  pipeline  optimized  for  MIPSGAL,  which  is  based  on 
the  Spitzer  Science  Center  (SSC)  pipeline  ('^S 1 6. 1 ).  It  includes 
the  standard  “read-2,”  row-droop,  and  droop  corrections,  dark 
subtraction,  flat  fielding,  and  a  modified  linearization.  The 
optimization  involves  better  artifact  mitigation  such  as  overlap, 
jailbar,  and  latent  corrections;  additional  details  are  described 
by  Mizuno  et  al.  (2008)  and  Carey  et  al.  (2009).  The  70  pm 
data  are  not  considered  here  and  thus  are  not  described.  For 
completeness,  we  note  that  BD+43  3710  was  outside  the  70  ^m 
coverage,  HBHA  4202-22  is  not  detected  in  the  post-bed  images, 
and  the  inner  ring  of  G79.29+0.46  is  apparent  in  the  post-bed 
images  but  the  new  outer  ring  is  not  (see  Section  3.3). 

The  Cygnus-X  survey  measurements  with  IRAC  consisted 
of  ~1?1  X  1°1  tiles  taken  in  high  dynamic  range  mode  with 
3  X  12  s  integrations.  The  IRAC  data  were  processed  using  the 
SSC  pipeline  SI 6.1  products  with  artifact  mitigation  performed 
for  jailbar,  pulldown,  muxbleed,  and  banding  effects  from  bright 
sources.  Mosaics  were  created  with  the  WCSmosaic  package 
which  performs  additional  corrections  detailed  by  Gutermuth 
et  al.  (2008).  Photometry  was  extracted  using  the  “Phot Vis” 
package  (Gutermuth  et  al.  2008)  for  the  IRAC  data  and  the 
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MIPSGAL  point-source  pipeline  (S.  Shenoy  et  al.  2010,  in 
preparation)  for  the  24  pm  data. 

BD+43  3710,  discovered  as  being  extended  in  the  MIPS 
data  taken  in  2007  December,  was  outside  the  IRAC  cov¬ 
erage  and  so  was  observed  in  a  separately  designed  Astro¬ 
nomical  Observation  Request  (AOR).  The  same  data  reduction 
process  was  applied  as  with  the  primary  IRAC  observations. 
HBHA  4202-22  was  at  the  edge  of  the  IRAC  coverage  region 
and  was  observed  only  with  the  3.6  and  5.8  pm  arrays,  as  the 
circumstellar  ring  was  not  discovered  in  time  to  adjust  the  AORs 
or  add  an  additional  field  as  was  done  for  BD+43  3710. 

We  also  observed  BD+43  3710  with  the  IRS,  targeting  the 
central  star  and  one  position  in  each  lobe  of  the  24  pm  bipolar 
nebula.  The  observations  used  all  four  low-resolution  modes  and 
the  data  (Pipeline  SI 8.0.2)  were  reduced  using  the  Spitzer  IRS 
Custom  Extraction  (SPICE)  software,  version  2. 1 .2,  provided 
by  the  SSC.  For  the  central  star  observations,  we  subtracted 
the  alternate  order  for  the  short- low  (SL)  data  and  the  alternate 
nod  for  the  long-low  (LL)  data  prior  to  performing  the  point- 
source  extraction.  After  the  spectra  were  extracted  with  SPICE, 
they  were  trimmed  to  get  rid  of  bad  data  at  the  extremes  of 
the  wavelength  ranges  and  averaged  together.  For  the  positions 
in  the  lobes,  we  used  the  full  slit  extraction,  which  includes 
corrections  to  the  calibration  due  to  the  aperture  size  and 
extended  nature  of  the  target.  There  were  no  compact  sources 
apparent  in  either  the  SL  or  LL  slit  for  the  lobe  observations; 
no  background  subtraction  was  performed.  The  luminous  blue 
variable  (LBV)  candidate  G79.29+0.46  was  observed  with 
IRS  in  high-resolution  mode  by  another  group  and  results  are 
presented  by  G.  Umana  et  al.  (2010,  in  preparation). 

Photometric  measurements  from  this  work  and  from  the 
literature  are  summarized  in  Table  3  for  each  of  the  three  stars. 
BVR  data  are  from  the  Naval  Observatory  Merged  Astrometric 
Dataset  (NOMAD;  Zacharias  et  al.  2004),  JHK  data  are  from 
the  Two  Micron  All-Sky  Survey  (2MASS;  Cutri  et  al.  2003), 
and  mid- infrared  data  are  from  our  Spitzer  measurements. 
The  uncertainties  for  the  Spitzer  data  are  the  extraction-based 
uncertainties  combined  with  the  absolute  uncertainties  for  IRAC 
(2%;  Hora  et  al.  2008a)  and  MIPS  (4%;  Engelbracht  et  al.  2007). 
Those  for  2MASS  are  from  Cutri  et  al.  (2003;  see  also  Skrutskie 
et  al.  2006);  the  BVR  uncertainties  are  estimated  at  4%  as  they 
are  not  listed  in  NOMAD. 

3.  RESULTS  AND  DISCUSSION 
3.1.  BD-^43  3710 

BD+43  3710,  located  in  the  extreme  northeast  comer  of  the 
Cygnus-X  survey  region,  is  a  little  known  carbon  star  candidate 
classified  as  spectral  type  R:  by  Nassau  &  Blanco  (1954).  It 
shows  a  remarkable  bipolar  structure  at  24  pm  that  is  completely 
absent  at  the  shorter  wavelengths  (it  is  also  outside  the  70 /xm 
MIPS  coverage).  Figure  1  shows  a  three-color  image  of  BD+43 
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Table  2 

Source  Parameters 


Name 

R.A.  (J2000) 

Decl.  (J2000) 

1 

b 

Spectral 

Type 

Object 

Type 

BD+43  3710 

20^45'"  34^3 

+43  3r2r28 

83.3755 

0.3400 

R: 

Star 

HBHA  4202-22 

20^17"^  8^:i2 

+41  o7'2r:o 

78.3203 

3.1536 

WR:,  AOI 

Em.  line  star 

G79 .29+0.46 

20*’31  "’42:28 

+40  ^IW.'! 

80.9005 

1.6435 

B[e]  I: 

Star 

Notes.  Object  types  as  reported  in  Simbad.  Spectral  types  are  from  Nassau  &  Blanco  (1954),  Kohoutek  &  Wehmeyer  (1999),  and  this 
work,  and  Voors  et  al.  (2()(X))  for  BD+43  3710,  HBHA  4202-22,  and  G79.29+0.46,  respectively,  although  see  the  text  for  discussion. 
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3710  combining  the  3.6  ^im,  S/xm,  and  24)Ltm  data,  while 
Figure  2  shows  the  individual  images  from  all  four  IRAC  bands 
plus  the  MIPS  lAfim  image.  The  lAfim  emission  is  roughly 
T.l  X  0.'9  or  ~  1 .34  X  0.45  pc,  assuming  the  Cygnus-X  distance 
of  1.7  kpc. 

The  central  star  is  in  the  Bonner  Dure hmu stem ng  ( Argelander 
1903),  Tycho  (H0g  et  al.  2000),  2MASS  (Cutri  et  al.  2003), 
and  Midcourse  Space  Experiment  (MSX;  Egan  et  al.  2003) 
point'source  catalogs.  However,  no  references  discussing  its 
physical  properties  could  be  found  in  the  literature.  The  visual 
magnitude  did  change  from  a  reported  9.2  mag  in  the  Bonner 
Durchmustemng  (Argelander  1903)  to  10.1  mag  in  the  early 
1990s  as  measured  by  Tycho.  Other  than  indicating  variability 
in  this  object,  which  is  not  an  unusual  feature  in  a  carbon 
star,  such  a  sparse  sampling  of  the  temporal  baseline  (albeit 
long)  does  not  allow  us  to  say  anything  meaningful  about  the 
pulsation  properties  of  the  star.  Figure  3  shows  the  spectral 
energy  distribution  (SED)  for  BD+43  3710,  as  well  as  the  other 
two  sources.  The  smooth  curves  are  Planck  curve  fits  to  the 
available  photometry.  A  single-temperature  graybody  docs  not 
well  represent  the  emission  from  any  of  these  sources,  although 
the  BVR  data  are  probably  suppressed  by  an  unknown  amount 
of  extinction. 

As  mentioned  above,  we  obtained  low-resolution  IRS  data 
on  BD+43  3710  and  one  position  in  each  lobe.  Figure  4  (left) 
shows  the  locations  where  the  spectra  were  taken.  The  right  side 


shows  the  resulting  spectmm  for  the  central  star.  The  sp>ectmm 
is  largely  featureless.  In  comparison  to  the  naked  stars  in  the 
ISOSWS  atlas  of  Sloan  et  al.  (2003),  it  is  most  similar  to  HD 
19557,  a  carbon  star  of  type  R  noted  by  Goebel  et  al.  (1983) 
as  having  particularly  weak  carbon  features  in  visible  and  near- 
IR  bands.  In  particular,  the  inflection  in  the  IRS  spectmm  at 
9  pm  is  well  matched  in  the  SWS  sp)ectmm  of  HD  19557, 
although  it  is  not  apparent  in  most  of  the  naked  carbon  stars 
in  the  SWS  atlas  (class  l.NC  of  Kraemer  et  al.  2002)  or  in 
the  other  dust-free  SWS  spectral  classes  (l.N,  l.NO,  l.NE, 
and  l.NM).  The  IRS  spectmm,  however,  does  not  show  the 
turnover  at  ~5.5  pm  expected  from  C3,  CO,  and  CN  absorption 
features  (Goebel  et  al.  1978)  typically  observed  in  carbon  stars 
(e.g.,  Aoki  et  al.  1998;  Zijlstra  et  al.  2006).  Several  hydrogen 
recombination  lines  are  detected  in  emission,  as  annotated  in  the 
figure.  The  strengths  of  the  Pf  a,  Hu  or,  and  H 1 8-7  do  not  change 
between  the  raw  data  and  the  background-subtracted  data,  unlike 
the  typical  nebular  lines  such  as  the  [Neii]  feature  at  12.8^tm, 
which  nearly  vanishes  after  background  subtraction.  Thus,  the 
hydrogen  recombination  lines  are  likely  associated  with  BD+43 
3710  and  its  nebula,  and  are  not  simply  in  the  foreground  or 
background  cloud.  Although  not  typical  for  carbon  stars,  Balmer 
lines  have  been  detected  in  emission  in  the  carbon  star  UV  Aur 
A  (Herbig  2009).  Thus,  we  find  that  the  carbon  star  spectral 
classification  for  BD+43  37 10  is  supported  by  the  IRS  spectmm 
but  not  definitively  confirmed. 
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Figure  2.  Individual  images  of  BD+43  37 10  in  the  IRAC  bands  (top  two  rows) 
and  in  MIPS  24  /^m  (bottom  image).  Grayscales  are  logarithmic. 

The  spectra  from  the  two  lobes,  shown  in  Figure  5,  are  quite 
similar  to  each  other,  with  a  cool  dust  spectrum  that  rises  with 
wavelength  past  the  end  of  the  IRS  band  and  a  few  broad  features 
and  low-excitation  fine-structure  lines.  The  SED  of  the  western 
lobe  rises  a  bit  more  steeply  with  wavelength  than  that  of  the 
eastern  lobe,  possibly  indicative  of  slightly  warmer  dust.  Neither 
spectrum  is  consistent  with  a  single-temperature  graybody,  but 
the  majority  of  the  dust  must  be  cooler  than  '^100  K  since  the 
SED  peaks  beyond  35 /im.  A  108  K  graybody  is  the  “best” 
fit  to  the  LL  data  for  both  lobes  and  is  shown  in  the  figure 
for  comparison.  As  can  be  seen  from  the  spectra,  unlike  in 
a  small  number  of  sources  where  the  circumstellar  emission 
in  the  MIPS  24  /im  band  has  been  found  to  arise  from  a  high 
excitation  [O  iv]  line  (e.g.,  Morris  et  al.  2004  (a  WN  star);  Morris 
et  al.  2006  (an  SNR  candidate);  Billot  et  al.  2009  (sources  of 
unknown  nature,  possibly  SNR  or  PNe  candidates);  N.  Flagey 
et  al.  2010,  in  preparation),  or  the  ~26-30/im  MgS  feature 
seen  in  some  post-AGB  stars  and  PNe  (e.g.,  Forrest  et  al. 
1981;  Goebel  &  Moseley  1985;  Hony  et  al.  2002;  Bemard- 
Salas  et  al.  2009),  the  nebular  emission  around  BD+43  3710 
arises  from  small  dust  grains.  They  are  probably  carbon-rich 
grains,  dominated  by  amorphous  carbon,  given  the  lack  of  either 
a  silicate  emission  feature  around  9.7  fim  or  silicon  carbide 
around  1 1  /um.  The  broad  features  present  in  the  spectra  at 
6.3,  7-8,  1 1.3,  and  16-1 8 /^m  typically  attributed  to  polycyclic 
aromatic  hydrocarbons  (PAHs),  are  not  likely  to  be  from  the 
BD+43  3710  nebula  since  there  is  no  corresponding  nebulosity 
in  the  IRAC  5.8  or  8.0  /um  image  that  can  be  readily  associated 
with  the  24 //m  lobes  (Figure  2).  The  PAHs  and  low-energy 


Figure  3.  SEDs  for  BD+43  3710  (top),  HBHA  4202-22  (middle),  and 
G79.29+0.46  (bottom).  Temperatures  for  the  best-fit  (minimal  x^)  graybody 
curve  (smooth  curves)  are  given  in  the  lower  left  of  each  panel.  Pluses  corre¬ 
spond  to  the  photometry  from  Table  3.  Uncertainties  are  smaller  than  the  symbol 
size. 

fine-structure  lines  seen  in  Figure  5  are  extremely  common  in 
active  star-forming  regions  such  as  Cygnus-X,  so  those  features 
likely  arise  from  the  fore/background  Cygnus-X  cloud. 

3.2.  HBHA  4202-22 

HBHA  4202-22  is  located  on  the  western  edge  of  the  survey 
region.  Kohoutek  &  Wehmeyer  (1999)  report  Hor  emission 
from  the  star,  citing  “Dolidze  (1971)  No.  3”  which  probably 
corresponds  to  Dolidze  (1971).  This  may  be  where  the  WR 
candidacy  given  in  SIMBAD  comes  from  although  we  could  not 
confirm  this  due  to  the  unavailability  of  the  Russian  circular,  and 
it  is  otherwise  unnoted  in  the  literature.  Figure  6  shows  a  three- 
color  image  from  the  3.6, 5.8,  and  24  /um  data,  and  the  individual 
images  are  in  Figure  7.  We  cannot  say  definitively  there  is  no 
emission  from  PAHs,  as  we  do  not  have  the  8  /um  data  which 
would  contain  the  (typically)  strongest  feature.  However,  there 
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Figure  4.  Left:  positions  (circles)  observed  with  IRS  in  BD+43  3710,  overlaid  on  the  24  data  (grayscale).  The  orientation  of  ihe  SL  and  LL  sliis  is  shown  in  the 
upper  left  comer.  Right:  IRS  5.25-21  nm  spectrum  of  BD+43  3710.  A  handful  of  hydrogen  recombination  lines  are  detected  in  emission  and  are  labeled  on  the  plot. 
The  feature  slightly  shortward  of  H  8-7  is  only  present  in  one  nod  and  so  is  probably  spurious. 


X  (pm)  X  (pm) 

Figure  5.  Left:  SL  IRS  spectra  from  the  two  lobes,  east  in  black  and  west  in  gray.  Right:  LL  IRS  spectra  frcrni  the  two  lobes,  east  in  black  and  west  in  gray.  The  labeled 
emission  features  likely  belong  to  the  fore /background  Cygnus-X  cloud,  not  BD+43  3710,  as  discussed  in  the  text.  A  108  K  graybody  is  the  "best”  fit  to  the  LL  data 
for  boih  lobes  and  is  shown  as  the  gray  dotted  line. 


is  little  extended  emission  in  the  5.8  /xm  data  that  would  have 
contained  the  6.2  /urn  feature  if  it  were  present  at  any  significant 
strength. 

Figure  8  shows  a  slice  through  the  ring  at  position  angle 
30°,  which  corresponds  to  the  narrowest  point  in  the  ring.  The 
sharper,  southwest  edge  of  the  ring  is  from  the  central  star, 
or  about  1.9  pc,  assuming  a  distance  of  6.5  kpc  (see  modeling 
discussion  below).  The  ring  is  more  extended  to  the  northeast, 
~1.'2,  or  about  2.3  pc,  and  the  edge  is  less  well  defined.  Two 
possible  explanations  come  to  mind.  First,  the  surrounding 
medium  could  be  less  dense  to  the  northeast,  allowing  the  shell 
to  expand  more  easily  in  that  direction.  However,  the  5.8  fim 
data  indicate  that  if  anything  there  is  more  material  to  the 
east-northeast  than  toward  the  southwest.  A  second  possibility 


is  that  the  star  is  moving  toward  the  southwest  and  the  material 
on  the  “front”  side  of  the  shell  is  l>etng  compressed.  Indeed,  the 
proper  motion  of  the  HBHA  4202-22  is  fia  =  —4.5  mas  yr“*, 
fis  =  —2.7  mas  yr“*  (Zacharias  etal.  2004),  which  corresponds 
to  a  position  angle  of  3 1  °,  nicely  consistent  with  the  compression 
direction.  The  distance  to  this  object  combined  with  the  proper 
motion  suggest  that  it  may  be  a  runaway  star. 

We  used  a  modified  version  of  the  radiative  transfer  code 
of  Egan  et  al.  (1988)  to  model  the  dust  shells  around  HBHA 
4202-22.  We  modeled  the  central  star  as  an  AO  supergiant  with 
r*  =  10,000  K,  L*  =  10^  Lq  at  a  distance  of  6.5  kpc,  i.e., 
well  behind  the  Cygnus-X  region.  While  these  parameters  are 
formal  inputs  to  the  model,  they,  too,  were  varied  in  order  to 
gel  the  best  fit  to  the  observed  data.  As  noted  above,  the  shell 


No.  6,2010 


CIRCUMSTELLAR  STRUCTURE  AROUND  EVOLVED  STARS  IN  THE  CYGNUS-X 


2325 


Figure  6.  Three-color  image  of  HBHA  4202-22.  Red:  24 /im;  green:  5.8  blue:  3.6 /im.  The  image  is  6^02  x  6!02  or  '^1 1.4  pc  at  a  distance  of  6.5  kpc  (see  the 
text). 


around  HBHA  4202-22  is  compressed  toward  the  southwest 
compared  to  the  northeast.  Therefore,  we  fit  two  models  to  the 
two  portions.  There  is  also  a  hint  of  a  ring  at  '^35"  in  the  5.8  fim 
data,  so  an  inner  shell  is  also  included  in  the  model.  Table  4 
gives  the  derived  parameters  using  amorphous  carbon,  with  a 
particle  size  of  0.135  /xm,  an  opacity  of  r8  035  fim  =1.5  x  10“^ 
(^0.55 /im  =  1.77  X  10“^),  and  grain  constants  from  Mathis  & 
Whiff en  (1988).  Figures  9  and  10  show  the  results  of  the  model 
compared  to  the  data  for  the  primary  shells  (compressed  and  not) 
at  24  fin\.  Models  using  silicate-rich  dust  were  also  considered, 
which  result  in  cooler  dust  grains.  The  match  to  the  data  was  not 
as  good,  though,  so  only  the  carbon-rich  results  are  presented  in 
Table  4. 

3.3.  079.29-^0.46 

Figure  1 1  shows  the  three-color  image  of  G79.29+0.46  from 
the  3.6,  8,  and  24 /xm  data.  It  was  first  found  to  have  a  ^^4' 
circumstellar  ring  in  the  radio  by  Higgs  et  al.  (1994),  who 
also  suggested  that  the  central  object  is  probably  an  LBV. 
Subsequently,  the  same  ring  was  detected  in  the  infrared  at  a 
number  of  wavelengths  from  8  to  60 /xm  with  IRAS  (Waters 
et  al.  1996),  the  Infrared  Space  Observatory  {ISO;  Wendker 
et  al.  1 998),  and  MSX  (Egan  et  al.  2002;  Clark  et  al.  2003).  This 
inner  ring  is  detected  with  all  four  IRAC  bands  as  well  as  the 
24 /xm  MIPS  band.  Figure  12  shows  the  individual  images  from 


Table  4 

Dust  Shell  Parameters:  HBHA  4202-22 


Primary  Shell 

Compressed 

Extended 

Inside  radius 

5.28  X  10**  cm 

5.28  X  10**  cm 

Thickness 

0.65  X  10**  cm 

3.4  X  10**  cm 

Dust  density  law 

r** 

Dust  temperature 

156.0-151.2  K 

156.0-137.4  K 

Dust  opacity 

X8.033  /tm  =  1.4  X  10 

ro.55/im  =  1.77  X  10~- 

Inner  shell 

Inside  radius 

5  X  \(y^  /?* 

3.70  X  10**  cm 

Thickness 

0.70  X  10**  cm 

Dust  density  law 

r  2 

Dust  temperature 

196.6-188.2  K 

Dust  opacity 

X8.035  /im  =  7.5  X  10^^ 

all  five  bands.  The  apparent  break  in  the  south  of  the  ring, 

particularly  in  the  IRAC  bands,  is  actually  due  to  an  infrared 
dark  cloud  (IRDC)  that  extends  south-southwest  in  the  figures. 
This  suggests  that  as  with  HBHA  4202-22,  G79.29+0.46  may 
be  more  distant  than  the  1.7  kpc  assumed  for  the  Cygnus-X 
complex.  Interestingly,  though,  the  CO  observations  of  Rizzo 
et  al.  (2008)  also  show  a  break  in  the  CO  emission  in  roughly 
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Figure  7.  Individual  images  for  HBHA  4202-22.  As  with  BD+43  3710,  there  is 
no  indication  of  extended  emission  associated  with  the  star  in  the  IRAC  bands 
(top  and  middle  panels),  in  contrast  to  the  circumstellar  shell  detected  at  24  /xm 
(bottom). 


the  same  position  as  the  IRDC.  Since  IRDCs  are  not  known 
to  absorb  molecular  line  emission  in  the  submillimeter,  this 
suggests  that  the  two  structures  might  actually  be  interacting 
in  some  fashion.  Additionally,  the  infrared  ring  appears  to 
be  somewhat  flattened  in  the  southeast,  also  consistent  with 
possible  interaction. 

Also  visible  in  the  24 /xm  data  in  Figures  11  and  12  is  a 
newly  discovered  outer  ring  ^1'  across.  This  large  ring  probably 
represents  an  earlier  episode  of  mass  loss  from  G79. 29+0.46 
compared  to  the  previously  known  inner  ring.  Additionally,  the 
24  /xm  emission  from  the  inner  ring  is  much  more  extended  than 
that  from  the  shorter  wavelengths.  This  kind  of  structure,  where 
the  IRAC  emission  ring  is  interior  to  a  larger  round  emission 
structure  at  24  /x  m,  is  occasionally  seen  in  the  MIPSG  AL  catalog 
of24/xm  rings  and  disks  (Mizunoet  al.  2010).  Usually,  the  8  /xm 
emission  appears  to  be  co-spatial  with  the  24 /xm  emission,  if 
present,  but  most  often  it  is  absent  entirely  (as  may  be  the  case 
with  HBHA  4202-22). 

As  with  HBHA  4202-22,  we  modeled  the  dust  shells  for 
G79. 29+0.46.  The  stellar  parameters  are  =  18,000  K, 
L*  =  4  X  10^  Lq,  =  4.567  x  10*^  cm,  with  a  distance 


Right  Ascension 


Figure  8.  Top:  24  /xm  ring  around  HBHA  4202-22.  Bottom:  slice  through  the 
30  position  angle  indicated  by  the  line  in  the  lop  figure.  This  angle,  which 
corresponds  to  the  narrowest  part  of  the  ring,  is  just  1  off  of  the  3 1  proper- 
motion  direction. 


HBHA  4202-22,  24  micron 


Figure  9.  Radial  profiles  for  HBHA  4202-22.  The  dots  show  the  24  data;  the 
dash-dotted  lines  show  a  horizontal  slice  through  the  data,  red  for  the  extended 
side,  blue  for  the  compressed  side.  The  solid  colored  lines  show  the  model 
results. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Figure  10.  Approximate  range  of  the  extended  and  compres.sed  shells.  Arcs 
showing  the  extent  of  the  dust  models  superposed  on  the  24  /im  image.  Their 
angular  extent  was  estimated  by  eye,  not  fit  by  the  model.  The  solid  horizontal 
line  indicates  the  radii  shown  with  the  red  and  blue  lines  in  Figure  9.  The 
dotted  line  with  the  arrow  indieates  the  direetion  of  the  proper  motion  of  HBHA 
4202-22. 


of  J  =  9.25  X  10^*  cm,  i.e.,  slightly  behind  Cygnus-X  at 
~3  kpc.  Table  5  lists  the  derived  parameters  for  the  inner  and 
outer  shells.  A  mix  of  amorphous  earbon  grains  and  PAHs  (Li  & 
Draine  2(X)1 )  was  used  in  the  model,  as  the  du.st  alone  could  not 
produce  the  flux  detected  in  the  IRAC  bands  and  the  strong  PAH 
emission  seen  in  the  IRS  observations  of  G.  Umana  et  al.  (2010, 
in  preparation).**^  The  amorphous  carbon  grains  dominate  the 
mass  and  opacity,  although  the  (much)  smaller  PAHs  are  more 
numerous.  While  the  inner  shell  is  modeled  as  a  single  shell 
with  a  density  law  of  the  24 /zm  flux  profile  (Figure  13) 
suggests  that  a  more  complex  model  with  a  different  density  law 
for  the  interior  of  the  shell  might  give  better  results. 

4.  SUMMARY 

We  have  detected  extended  24 /zm  emission  around  three 
evolved  stars  in  the  Cygnus-X  Legacy  survey  region.  A  bipolar 
nebula  was  found  around  the  earbon  star  BD-h43  3710.  IRS 
spectra  show  that  the  emission  is  due  to  small  dust  grains, 
rather  than  [Oiv]  line  emission  or  the  broad  MgS  feature. 
An  asymmetric  ring  of  24  /zm  emission  was  observed  around 
the  Wolf-Rayet  candidate  HBHA  4202-22.  The  '‘compressed** 


Morris  et  al.  (2008)  present  an  infrared  spectrum  of  G79.29+0.46  that 
included  IRS  dala  as  well  as  ISOSWS  and  near-infrared  spectra.  Presumably 
ihe  IRS  spectrum  is  from  the  central  star  as  it  does  not  show  the  strong  PAHs 
lhal  are  detected  by  G.  Umana  et  al.  (2010.  in  preparation)  toward  the  shell. 


Figure  11.  Three-color  image  of  G79.29+0.46.  Red:  24  ^m;  green:  S/zm;  blue:  3.6 /xm.  The  image  is  15'  x  15'  or  ''-8.75  pc  at  a  distance  of  2  kpc,  or  13. 13  pc  at  a 
distance  of  3  kpc  (see  the  text). 
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Figure  12.  Individual  images  of  G79.29+0.46  in  the  IRAC  bands  (top  two  rows) 
and  at  24  fjim  (bottom).  Note  that  the  gap  in  the  south  of  the  inner  ring  is  caused 
by  foreground  absorption  associated  with  the  IRDC  stretching  across  the  bottom 
of  the  images.  Grayscales  are  logarithmic. 


portion  of  the  ring  is  in  the  direction  of  the  projected  velocity 
vector  derived  from  the  observed  proper  motion  of  the  star. 
Radiative  transfer  models  imply  that  the  central  star  is  an  A-type 
supergiant  at  a  distance  of  '^6.5  kpc  and  the  ring  is  composed 
of  amorphous  carbon  grains.  A  large,  '^1'  diameter  ring  was 
detected  at  24 /xm  around  the  LBV  candidate  G79.29+0.46, 
a  new  addition  to  the  well-known  inner  ring.  The  inner  ring 
was  detected  at  all  four  IRAC  bands  as  well  as  at  24  /xm.  The 
shells  were  modeled  with  a  mixture  of  amorphous  carbon  grains 
and  hot  PAHs.  An  IRDC  blocks  part  of  the  ring  emission  at 
the  shorter  wavelengths,  which  along  with  the  model  results 
indicates  that  the  LBV  candidate  is  at  least  on  the  far  side  of  the 
Cygnus-X  region,  if  not  further. 

While  somewhat  disparate  in  their  properties,  these  three 
objects  demonstrate  the  potential  of  the  Spitzer  data  archives. 
The  thrust  of  the  Cygnus-X  Legacy  project  was  to  explore  the 
richest  high-mass  star  formation  region  within  2  kpc.  Similarly, 
the  larger  Spitzer  Legacy  programs  that  mapped  the  inner 
Galactic  plane  with  IRAC  (GLIMPSE,  the  Galactic  Legacy 
Infrared  Mid-Plane  Survey  Extraordinaire;  Benjamin  et  al. 
2003)  and  MIPS  (MIPSGAL;  Carey  et  al.  2009)  also  emphasized 
star  formation  studies  in  their  proposals  and  overview  papers. 
In  each  case,  though,  a  wealth  of  information  on  circumstellar 
structures  around  evolved  stars  also  has  been  obtained.  For 
example,  hundreds  of  small  24 /xm  emission  rings  have  been 
found  in  the  MIPSGAL  survey  data  (Mizuno  et  al.  2010),  many 
of  which  are  similar  to  the  structures  around  HBHA  4202-22 


LBV  G79,  24  pm 
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1000 
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Figure  13.  Radial  profiles  of  the  rings  around  G79.29+0.46.  The  red  line  is 
a  slice  through  the  data  westward  of  the  central  star.  The  blue  line  shows  the 
model  fit. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


Table  5 

Dust  Shell  Parameters:  G79.29+0.46 


Parameters 

Inner  Shell 

Ouier  Shell 

Inside  radius 

3.973  X  10**  cm 

8.99  X  10**  cm 

8.7  X  10^/?, 

LI  dinner 

Thickness 

4.2  X  10**  cm 

0.5  X  10**  cm 

Dust  densiiy  law 

^-3.5 

r-2 

Dust  temperature 

Am.  carbon 

107.6-89.9  K 

87.75-86.57  K 

PAHs* 

1500-354  K 

96.1-94.4  K 

Dust  components^ 

^^irnxw. 

0.45  pm 

0.40 /im 

^^PAHs 

5.0  X  10*"^  pm 

5.0  X  10  pm 

2% 

3% 

^^PAHs 

98% 

97% 

Dust  Opacity 

^24  //m 

1.324  X  10-^ 

9.81  X  10^^ 

fs.o  pm 

2.70  X  10-^ 

2.00  X  10 

ro.55  pm 

1.27  X  10-2 

9.39  X  10 

Notes. 

*  PAH  temperature  is  non-equilibrium  for  the  inner  shell  and  declines  as  r 
The  outer  shell  is  an  equilibrium  temperature. 

^  a  =  grain  size;  N  =  number  density  within  the  shell. 


and  G79.29+0.46  (the  main  difference  being  that  most  of  the 
MIPSGAL  rings  have  no  detected  24 /xm  central  source).  As 
with  our  two  ring  sources,  those  objects  that  have  previous 
identifications  are  often  identified  as  emission  line  stars.  It  is 
also  interesting  to  note  that  we  detect  extended  emission  in  the 
IRAC  bands  only  toward  the  inner  ring  of  G79.29+0.46,  and 
the  MIPSGAL  rings  have  IRAC  counterparts  less  than  15%  of 
the  time.  Data  from  follow-up  IRS  observations  are  currently 
being  analyzed  to  determine  what  mechanism  is  responsible  for 
the  24 /xm  emission  in  a  selection  of  the  MIPSGAL  rings  (N. 
Flagey  et  al.  2010,  in  preparation) — amorphous  carbon,  as  in 
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BD+43  3710;  [Oiv]  as  in  the  Morris  et  al.  and  Billot  et  al. 
SNRs;  the  as-yet-unidentified  21  /zm  feature  seen  in  post-AGB 
stars;  or  some  other  mechanism.  Results  from  these,  as  well  as 
other  studies  underway,  will  enable  us  to  better  place  our  three 
objects  in  the  broader  context  of  evolved  stars  with  circumstellar 
envelopes  in  the  Galaxy. 

Note  added  in  proof.  We  note  the  following  papers  were 
published  during  the  review  process  which  may  be  of  interest 
to  the  reader,  as  they  discuss  additional  aspects  of  G79.29+0.46 
and  HBHA  4202-22,  respectively:  Jimenez- Esteban  et  al.  2010, 
ApJ,  713,  429  and  Gvaramadze  et  al.  2010,  MNRAS,  in  press; 
and  Gvaramadze  et  al.  2009,  MNRAS,  400.  524. 

This  work  is  based  in  part  on  observations  made  with  the 
Spitzer  Space  Telescope,  which  is  operated  by  JPL/Caltech 
under  NASA  contract  1407.  Support  for  this  work  was  provided 
in  part  by  NASA.  This  research  has  made  use  of  NASA’s 
Astrophysics  Data  System  Bibliographic  Services,  data  products 
from  2MASS  which  is  a  Joint  project  of  the  University  of 
Massachusetts  and  IPAC/Caltech  funded  by  NASA  and  the 
NSF,  and  the  Simbad  database  operated  at  CDS,  Strasbourg, 
France.  We  thank  the  anonymous  referee  for  useful  suggestions 
that  improved  the  paper. 

Facilities:  SpitzerifRNC,  IRS,  MIPS) 
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